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ABSTRACT 

We study steady, radial gas outflows from galaxies in an effort to understand the way tenuous and 
hot gas is transported to large distances away from galaxies. In particular, we obtain solutions for outflow 
problems and study the outflow topology, the effects of the galaxy potential (and thus the mass), the size of 
outflow regions, the efficiency of radiative cooling, and the fate of the cooled gas. Under general power-law 
forms for the cooling function and the gravitational field of galaxies, we show that the outfiow solutions are 
determined by the two-parameter initial conditions. In an analogy with stellar wind or accretion problems, 
we demonstrate that there exist no transonic fiows, but either subsonic or supersonic fiows are obtainable. 
Solutions of the supersonic outfiows are studied in detail as they are most likely to carry gas to large distances 
away from galaxies. We find that, if gravity is weak, the outfiow is characterized by the ratio of the radiative 
cooling time to the fiow time, tc/tj. If initially tc/^/ ~ 1 the gas cools as soon as it leaves the galaxy, whereas 
if tc/tf ^ 1 the gas first cools adiabatically and then radiatively. However, if initially tc/tf 1 radiation 
cannot become important; in this case, depending on whether the initial gas temperature is above the escape 
temperature, the outfiow results in either a galactic wind or a hot corona. The importance of the galactic 
gravitational field is characterized by the fractional energy lost radiatively within the fiow time in outfiows 
with velocity equal to the circular velocity of the galaxy; if the fraction is small, gravity is relatively strong, 
and it stops the outfiow before the gas has a chance to cool radiatively, resulting in a hot corona. In case the 
gas does cool radiatively, the cooled gas is most likely to form clouds via various instabilities. The clouds 
coast farther away from the galaxy because of the finite kinetic energy they inherit. Depending on the initial 
energies, the clouds can either leave the galaxy or fall back ballistically. We apply the calculations to galaxies 
in general. For reasonable ranges of the parameters, we find that the hot gas in dwarf galaxies can either fiow 
out as galactic winds, or cool radiatively to form clouds. In the latter case, the clouds escape the galaxies. 
In contrast, massive galaxies like our own tend to confine the gas. The gas released into the halo can either 
cool radiatively or result in a galactic corona. We present the surface brightness in various x-ray energy 
bands for some representative cooling outfiows from dwarf and normal galaxies, and we find that the extent 
of the resultant X-ray emission is generally much smaller than the region of the outfiows. In particular, the 
mean temperature averaged over entire outfiow regions is shown to be a factor of 2-6 smaller than the base 
temperature. We also estimate the mean surface brightness of the O VI emission lines, and the predicted 
surface brightness is within the reach of current UV experiments. We briefiy discuss the implications of cold 
clouds at large distances from dwarf galaxies for recent observations of the QSO absorption line systems. 

1. INTRODUCTION 

Studies of gas outfiows from galaxies are of considerable importance to understanding various galaxy 
phenomena. For example, the ejection of enriched gas dilutes the abundance of the interstellar gas in galaxies, 
thus affecting the subsequent rate of radiative cooling and star formation (e.g., Larson 1974; Wang & Silk 
1993). Outfiows can transport momentum and energy out of galaxies, therefore changing the dynamics of 
the galaxies (e.g., Dekel & Silk 1986). The heavy elements and kinetic energy injected from galaxies through 
outfiows may have caused heating and enrichment of the intracluster medium and intergalactic medium in 
the early universe (e.g.. White 1989). 

Of more recent interest is the increasing realization that the gas found through QSO absorption lines 
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such as the Mg II lines is associated with galaxies at extremely large distances. Indeed, attempts to identify 
spectroscopically the underlying galaxies at intermediate redshifts have been successful, indicating that the 
Mg II absorbing gas is at a typical distance of about 50 kpc from the galactic centers (e.g., Bergeron & 
Boisse 1991; Yanny et al. 1990; Bechtold & Ellingson 1992; Steidel 1993; Steidel & Dickinson 1995). More 
recently, studies of the Lya clouds in relation to nearby galaxies suggest that they may be also associated 
with galaxies, possibly at distances even larger than the heavy element absorption line gas (Morris et al. 
1993; Lanzetta et al. 1995). At high redshifts, observations of pairs of quasars show that the Lya clouds 
have an average radius of about 100 kpc (Bechtold et al. 1994; Dinshaw et al. 1995). 

A possible explanation is that the absorbing gas results from gas outflows from the main body of the 
galaxies where most of stars and gas reside (Bregman 1981). Mathews & Baker (1971) first considered 
outward-fiowing galactic winds from elliptical galaxies as a result of the gas ejection from stars. However, 
here we are mostly concerned with late-type (star-forming) galaxies, both because they are generally more 
numerous and because active star formation is inferred in the absorbing galaxies (e.g., Steidel 1993). In star- 
forming galaxies, the interstellar gas is shock-heated by stellar winds and more importantly, by supernova 
remnants to an extremely high temperature (e.g.. Cox & Smith 1974). That the resulting hot gas may 
fiow out of the galaxy was first proposed by Shapiro & Field (1976), and the possibility that the gas cools 
radiatively in the halo and subsequently falls back to the disk ballistically is called a "galactic fountain" . 

Subsequent hydrodynamic calculations of the galactic fountains have focused mostly on gas fiows in 
galactic potentials similar to our Galaxy (Bregman 1980, 1981; Habe & Ikeuchi 1980; Houck & Bregman 
1990; Li & Ikeuchi 1992; also see Cox 1981). However, simple arguments show that in galaxies like our 
own the heating rate required for the gas to reach a distance from the galactic center comparable to that 
observed in the metal line systems is implausibly large (Wang 1993). This is not surprising because the 
deep gravitational potential of galaxies requires extremely energetic gas to move out of the potential. This 
perhaps also explains the scarceness of the large X-ray halos detected for normal spiral galaxies like our 
own by ROSAT. However, for less massive (dwarf) galaxies one expects that gas can reach large distances 
because of the galaxies' shallow potential wells (e.g.. Silk & Norman 1979; Fransson & Epstein 1982; Dekel 
& Silk 1986, Vader 1987). 

Despite the increasing realization of the importance of gas outfiows from galaxies, many questions 
concerning the physical nature of outfiows remain unanswered. In particular, the dynamics of the outfiows, 
such as the topology, the condition for radiative cooling to become important, the effects of gravitational 
field (thus the galaxy mass), and the extent of outfiows, has not been studied systematically. As the first in 
a series to study gas outfiows from galaxies, the present paper is aimed at obtaining exact solutions for the 
outfiow problems and answering the above questions in a systematic way. 

The temperature and volume filling factor of the hot phase in the interstellar medium (ISM) resulting 
from supernova heating depends on the supernova rate per unit volume and the density of the ISM, which 
are not precisely known even for our Galaxy. Thus we shall assume that the outfiow originates from a 
fixed initial (base) radius at an initial temperature; the heating of the gas in the ISM is not considered 
here. Instead we consider ranges of the initial temperature and the initial radius as representing a variety of 
galaxies. If the porosity of the hot phase in the ISM is close to unity, the base radius is roughly the edge of 
the stellar distribution where the heating by supernovae decreases substantially. If the interstellar volume is 
dominated by cold clouds, the base radius should then be roughly the radius where the hot gas breaks out of 
the cold medium (such as in the chimney model; Norman & Ikeuchi 1989). We ignore thermal conduction, 
which may be unimportant in the large-scale dynamics of the upward fiow since the heat fiux conducted into 
the halo is usually smaller than the energy radiated (Houck & Bregman 1990). Any possible magnetic field 
(ignored here) will further impede thermal conduction. 

We shall consider only steady outfiows. Typically, the mass of the hot phase in the ISM is negligible 
compared with the total mass of the ISM, so a sufficient gas reservoir is generally available. The the steady 
fiow calculations can be applied to galaxies where the duration of energy injection from supernovae into 
the ISM is sufficiently long, so the temperature at the base of the halo may be kept fixed during the fiow, 
and a steady state can be achieved. For a velocity of 150 km/s (corresponding to the sound speed in a gas 
of 10® K), the gas advances 15 kpc within 10® years. A burst of star formation, even it is a (5-function in 
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time, will result in the gas heating with a duration of at least a few times 10 years for the standard initial 
mass function because of the finite lifetime of the massive stars that eventually explode as supernovae (e.g., 
Leitherer & Heckman 1995). Thus the steady fiow solutions can be applied to dwarf galaxies where the 
outfiow region is typically less than a few kpc (see below), regardless of whether star formation is in a burst 
or in a continuous phase. They are also applicable to normal disk galaxies like our own, where star formation 
has been continuous and has varied little for the last few billion years (e.g., Sandage 1986; Wang & Silk 
1994). Finally, we limit our calculations to radial (spherically symmetric) fiows, which may be a reasonable 
approximation if the base of the fiow is much smaller than the region of the outfiow. 

In §2 we first show that for power-law and logarithmic potentials where no additional scales are in- 
troduced, the outfiow equations can be reduced to dimensionless form by introducing two natural scales 
in the problem (§2a). We then demonstrate that the outfiow solutions cannot be transonic, in an analogy 
with stellar wind or accretion problems (§2&; Appendix), and that they can be obtained in terms of the 
two-parameter initial conditions. We then give examples of the supersonic solutions (§2c), and discuss the 
general properties of the outfiows (§2(i). In §3, we apply our calculations to dwarf and normal galaxies. 
Finally, we present a discussion and summary in §4. 



2. COOLING OUTFLOWS 



2a). Formulation of the Problem 
The equations of mass, momentum, and energy conservation in a steady, spherically symmetric, fiow 

1 d ir'^pv) 

-t: ^/ ^ = (2.1) 



dr p dr dr 

'^.ln(^\=-C, (2.3) 



7 — 1 (ir 

where p, v and p are, respectively, the density, velocity, and pressure of the gas, 7 (we take 7 = 5/3 unless 
otherwise noted) is the adiabatic index, (f) is gravitational potential (self-gravity of the gas is not included), 
and C is the cooling rate per unit volume. We ignore heating from possible intergalactic UV background 
radiation since it is in general not important for the hot gas considered here. The local adiabatic sound 
speed is 

/ \ 1/2 

..(!£) ,M 

As P is proportional to the gas temperature, Tg, the cooling rate per unit volume can then be written as 

C = Ap^c^i. (2.5) 



For example, for a gas of near-cosmic abundance of heavy elements, cooling is dominated by bremsstrahlung 
at Tg ^ 5 X lO^A' so q = 1/2. At the lower temperatures, cooling is dominated by the line emissions, 
and q ~ —1/2 to —1 (Dalgarno & McCray 1972; Raymond, Cox, & Smith 1976; Gaetz & Salpeter 1983; 
Sutherland & Dopita 1993). Using equation (2.1) to give 

Airr^pv = A = const., (2.1a) 

where A is the (constant) mass-loss rate, and substituting p in equations (2.2) and (2.3), we obtain 

dv^ _ (1 - -/)Ac^iX - Sirrc^v^ + iTTr^v^(d(f)/dr) 
dr 2nr'-^ (c^ — v'-^) 
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dc^ _ (1 - 7) [(c^ - jv^)AP''X - Snrv'^P + Anr^ P {d<j, / dr)] 
dr Airr^v^ (c^ — v^) 

It is thus clear from the above equations that the sonic point, where = P , is the critical point where 
the numerator vanishes, and the derivatives of v and c become infinite unless the numerators in equations 
(2.6)-(2.7) also vanish. If we use the Mach number, M = v/c, as the variable (e.g., Bondi 1952), then the 
critical point is at = 1. Notice that the above equations are invariant under the transformation v —v 
(or M —M), so the equations describe both gas outfiow and accretion. 

We now consider some specific forms of galactic gravitational potential, which is related to the circular 
(or rotation) velocity of the galaxy, v^r, by 



For simplicity, we assume a rotation curve 



r^-^r- (2.8) 



^cV = ^- (2.9) 



For example, for Keplerian orbit, n = 1 and B = GMq, where G is the gravitational constant, and Mq is the 
total mass interior to the orbit. For galaxies dominated by dark halos with a fiat rotation curve (isothermal 
spheres), we have n = and B = v^^^ = constant. 

Given the gravitational potential in equation (2.9), there exist only two dimensional quantities in the 
present fiow problem specified by equations (2.6)-(2.7): a characteristic length, ro, and a characteristic 
velocity, v^} It is easy to show, from dimensional arguments using A, B, and A, that these characteristic 
scales can be defined as the following: 



AX 



l/(l-2n+nq) 



'°=i^) (2-10) 



B 



l/2(l-2n+nq) 



^^=[a^) • (2.11) 

The characteristic radius ro and velocity vq above contain information concerning the radiative cooling and 
the gravitational field of the galaxy. This is most obvious if we consider isothermal spheres (n = 0): equations 
(2.10) and (2.11) now read 

ro = (2.10a) 

V ■ 

cir 

Vq = Vcir- (2.11a) 

Thus in this case, ro specifies the rate of the radiative cooling in a fiow with velocity equal to the rotation 
velocity of the galaxy, and vq is the rotation velocity itself, determining the strength of the gravitational 
field. Note that the mass-fiow rate enters here because it determines the density of the gas and therefore the 
cooling rate. 

We can now define dimensionless variables in place of the radius and sound speed (or temperature) : 

x=—; (2.12) 

ro 

w=—. (2.13) 

The dimensionless radius, x, is a parameter measuring the strength of the gravitational field of the galaxy 
relative to the radiative cooling rate (just as ro is). In fact, if n = we have 

2 

(2.12a) 



4:TT(r/Vcir)C 



^ The gas density can be found from mass conservation once solutions of (2.6)-(2.7) are obtained. The 
density scales as po = A/(47rroi'o), and the pressure scales as po = t;oA/(477rro). 
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That is, 1/ X is the fraction of the energy lost to radiation during the flow time r jv^r for a flow with velocity 
equal to the rotation velocity of the galaxy. Therefore, if a; ^ 1 one expects that gravity becomes important 
in the flow. The dimensionless temperature, w^, is a parameter characterizing the gas thermal energy relative 
to the potential energy. Thus for 1 we expect that the gas escapes the galaxy. 

Using the Mach number instead of the velocity, and substituting equations (2.12) and (2.13) into equa- 
tions (2.6) and (2.7), we obtain the dimensionless equations for the Mach number and the dimensionless 
temperature: 



dM'^ _ (7 - 1) (1 + 7M2) + 87ra;i+"w4^2 [2 + (7 - 1) M'^\ - 47r (7 + 1) xw'^ M'^ 

dx ~ 47ra;2+nu,4(7\^2 _ X) 

du? _ (7 - 1) {-x'^w'^'i {-fM"^ - 1) - 87ra;i+"w4^4 _^ ^t^xw'^m'^] 
~d^ ~ 47ra;2+nu,2 7\^2 {^^2 _ 1) • 



(2.14) 



(2.15) 



The first, second, and third terms in the numerators of equation (2.14) are due to radiative cooling, the 
gas pressure, and gravitation, respectively. Similarly, the first, second, and third terms in the numerator 
of equation (2.15) are contributions separately from radiative cooling, adiabatic expansion, and gravitation. 
The reduction to dimensionless equations above greatly simplifies our investigation of the fiow solutions in 
relation to various galaxy properties (such as the mass-loss rate, the temperature at the base of the halo, the 
depth of the galactic potential, and so on). The complexity of taking into account varying galaxy parameters 
is reduced to a set of initial condition problems. 

Two crucial timescales in the present problem are the fiow time and the cooling time. We define the 
fiow time as 

t; = ^, (2.16) 

and the cooling time as 

The ratio of the cooling time to the fiow time in terms of our dimensionless variables is: 

tr 47r 



if 7(7 - 1) 



-xw 



(2.18) 



In fact the above ratio is roughly the ratio of the contributions from adiabatic expansion and radiative 
cooling in supersonic fiows described by equations (2.14) and (2.15) (i.e., the second and first terms in the 
numerators). Iftc/tf ^ 1, the fiow is adiabatic, and iftc/tf ^ 1, radiative cooling is important, and the gas 
cools rapidly to low temperatures. 

2b). Supersonic Flow 



In the absence of radiative cooling (i.e., q —00), and for the point-mass potential (i.e., n = 1), 
equations (2.14) and (2.15) reduce to the classic Bondi problem (Bondi 1952); Bondi solved a spherical 
accretion problem for fixed adiabatic indices ranging from 7=1 (isothermal) to 7 = 5/3 (adiabatic). The 
corresponding wind solutions have been discussed by Parker (1958) in the context of the hydrodynamic 
model of the solar wind. From these classic studies, it is well understood that solutions for transonic fiows 
exist only for 7 < 5/3. There are no solutions proceeding smoothly from a subsonic to a supersonic regime 
owing to the stiff pressure response of an adiabatic gas; the sonic point progressively moves toward the origin 
as 7 — 5/3. Solutions do exist, however, for either subsonic or supersonic fiows for 7 > 5/3. The branch of 
the subsonic solutions turns out to require a finite confining pressure at infinity (Parker 1958). 

For general potentials with arbitrary values of n, the topology of the solutions with constant 7 differs 
from the n = 1 case. In the Appendix, we present fiow solutions for spherical accretion and wind problems at 
a fixed 7 (but without the cooling term included) in general potentials given by equation (2.9) with arbitrary 
values of n. Our main results there can be summarized as the following: (1) The range of 7 within which a 
transonic fiow is possible is progressively reduced as n decreases below unity. In particular, for isothermal 
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spheres (n = 0), a transonic flow can occur only if 7 < 1. Obviously, as the strength of the gravitational 
held increases toward small values of n, the gas flnds it more and more difficult to overcome the sonic barrier 
unless heating is extremely effective. (2) For a given value of n, there exists a critical value of 7 (e.g., 
7 = 5/3, 1 for n = 1, 0, respectively), above which flow can only be either subsonic or supersonic. If the 
total energy E > 0, the gas can reach to inflnity, with the Mach number increasing (decreasing) toward 
large radii for supersonic (subsonic) flows. However, in subsonic flows the density and temperature approach 
asymptotically constant values, implying a flnite conflning pressure at inflnity, a generalization of the result 
found by Parker (1958) for the n = 1 case. (3) For supersonic or subsonic flows in the above case, if the total 
energy E < Q, there is a maximum radius beyond which no flow solutions exist. For n ^ 0, the flow starts 
either as subsonic and advances to the sonic point or as supersonic and abates to the sonic point, as the flow 
approaches the flnal radius where it is stopped by gravity. For isothermals spheres (n=0), however, there 
are two sonic points. For outflow problems, the gas has to start at the inner sonic point in both supersonic 
and subsonic flows, and it approaches the outer sonic point as the flow reaches the flnal radius. 

With radiative cooling (q = flnite) (and 7 = 5/3), pressure response to density variations is even steeper 
than in the adiabatic flow. This effectively increases the adiabatic index, and from the Appendix we see that 
there exists no transonic flow solution given the potential deflned in equation (2.9). This indeed is borne 
out by our numerical analyses of equations (2.6)-(2.7). Solutions do exist, however, for either supersonic or 
subsonic flows. A potential problem with the subsonic flows is the outer boundary conditions mentioned 
above. Furthermore, because the velocities are relatively small in subsonic flows, the flows are close to 
hydrostatic equilibrium, as can be inferred directly from equation (2.2). Thus gravity is always important 
in subsonic flows. The entropy of the gas decreases toward large radii as a result of radiative cooling. 
Consequently, the subsonic flow is subject to convective and thermal instabilities (Balbus & Soker 1989). 
Hydrodynamic simulations of subsonic cooling outflows indeed bear this out (e.g., Houck & Bregman 1990). 
The supersonic flows, however, are driven by thermal pressure, and gravity is relatively unimportant for 
most of the flow. The gas becomes unstable only when it cools radiatively and gravity becomes important 
at the end of the flow. Various hydrodynamic simulations of supersonic cooling outflows show this is indeed 
the case (Bregman 1980; Habe & Ikeuchi 1980; Li & Ikeuchi 1992). 

Aside from instabilities and possible outer boundary problems, subsonic flows have smaller velocities 
than their supersonic counterparts, and the resultant flow densities are higher for given mass-loss rate and 
initial gas temperature. Thus the radiative loss is more efflcient and the resulting halo size is smaller than 
in the supersonic flows. As we are mostly interested in the largest extent of the galactic halos resulting from 
outflows, below we shall concentrate on the supersonic flows. Furthermore, we shall consider isothermals 
spheres (n = 0) for speciflc numerical examples since many galactic potentials can be well approximated by 
a flat rotation curve. 

The initial Mach number of the flow is taken as unity at the base of the flow, a natural choice for the 
initial condition in supersonic (or subsonic) flows. Indeed, in a calculation of winds from starburst galaxies. 
Chevalier & Clegg (1985) obtained an analytical solution for the wind driven from a region of uniform mass 
and energy deposition; both gravitational held and radiative cooling were ignored but a heating term was 
added to the flow equations (cf. eqs.[2.1]-[2.3]). These conditions may be applicable to the region of the 
gas heating (the base); here the heating of the gas is more important than radiative cooling, a prerequisite 
for the existence of the hot phase in the ISM in the flrst place. Furthermore, the gravitational potential 
within the base should not vary largely because of the relatively small size of the heating region, and also 
because the potential is dominated by stellar populations rather than dark matter (so the variation of the 
potential is not as steep). Chevalier & Clegg found that a smooth transition from subsonic flow at the center 
to supersonic flow at large radius r requires that M = I at r = R, where R is the radius beyond which the 
heating ceases. The conditions imposed in the present work (i.e., inclusion of gravitational held and radiative 
cooling) are appropriate for the gas beyond the heating region, and thus our solutions can be considered as 
an extension of their inner (r < R) solutions; the transition is at r = i? with M = I. 

We thus see that solutions of equations (2.14) and (2.15) are determined by two initial conditions Xi 
and wf at the base radius. These two parameters reflect the initial conditions of the gas, namely the relative 
importance of the gravitational potential to the radiative cooling, and of the internal thermal energy to the 
potential energy. Therefore, variations in Xi and wf represent the different conditions of galaxies such as the 
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region of the gas heating, the gas temperature and the mass of the galaxies. 

2c). Examples of Outflow Solutions 

As examples we show solutions for supersonic outflows with several values of the initial dimension- 
less temperature w'f and the initial dimensionless radius Xi in figures la-lc. The adopted values are [xi, 
w?) = (10~'*, 25) (figure la), (0.6, 1) (figure lb) and (10, 1) (figure Ic). Because equations (2.14) and (2.15) 
diverge at M'^ = 1, in practice we start the outward integration at M'^ = 1 -|- e, with e ^ 1 (typically 
10~^) being sufficiently small so that backward integration to M'^ = 1 does not in practice change any initial 
conditions imposed. The cooling exponent is assumed to be g = —0.6, corresponding to a gas of nearly 
cosmic abundance at a temperature 10® ^ Tg ^ 4 x 10^ K (Dalgarno & McCray 1972; Raymond, Cox, & 
Smith 1976; Gaetz & Salpeter 1983; Sutherland & Dopita 1993). Gravitational field of equations (2.8) and 
(2.9) with n = is adopted. 




The two sets of (xi, w'f) in figure la and lb have a similar ratio of the cooling to the fiow time, tc/tf , 
despite their large difference in the initial conditions. The values of [xi, wf) in figure Ic on the other hand 
lead to a much larger tc/tf even though the initial dimensionless temperature wf is the same as in figure 
lb; for given wf the relative importance of radiative cooling is reduced by increasing Xi (e.g., by lowering 
the mass-loss rate). As can be seen in the figures, for the parameters considered here radiative cooling is 
initially unimportant in all the cases and the gas first cools almost entirely owing to adiabatic expansion. In 
figures la and lb, as gas streams out toward large radii, radiative cooling eventually becomes important and 
the gas cools rapidly. In figure Ic, however, the time for radiative cooling is so long compared with the fiow 
time that the fiow remains adiabatic; here, before the gas has a chance to cool radiatively, the gas is pulled 
back by the strong gravitational field (characterized by Xi > 1), a situation already seen in the analytical 
solutions in the Appendix (cf. figure A3b). As the gas slows down and the density increases, the gas is 
in fact heated up through adiabatic compression. Because the fiow time is much shorter than the cooling 
time in figure Ic, the radiative cooling term can be ignored. We then expect that the gas in the halo should 
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reach a quasi-hydrostatic equilibrium over a few crossing times, resulting in a hot corona. Given enough 
time (many crossing times), the gas will cool gradually, with the cooled gas falling back to the galaxy. The 
amount of mass cooled is simply the lifetime times the cooling rate and should be much smaller than the 
total mass in the hot corona as long as the cooling time remains long. 




The primary driving force of the flow is thermal pressure; the gas is first accelerated by the gas pressure, 
which is much larger than gravity, and the thermal energy is converted to the kinetic energy. A unit mass 
of fiuid carries with it an amount of internal energy Ei = h + /2 (excluding the potential energy), where 
h is the enthalpy. Thus initially the gas has an energy per unit mass of 

E.= ^J^'^_\^ c^ = 2wfvl (2.19) 

where c,- is the initial adiabatic sound speed. Part of the this energy has to be used to overcome gravitation, 
so the maximum velocity is less than \/2Ei = 2wiVo. The gas cools as it expands, and the pressure drops. 
Eventually gravitation becomes dominant, and pulls the gas back. The actual peak of the velocity is reached 
roughly when the pressure force is equal to gravity; from equations (2.6) and (2.7) we find this occurs when 
~ 1/2 for a purely adiabatic fiow. The energy loss to radiation can be estimated from the energy difference 
between the initial and final states of the fiow. Iftc/tf is initially small, radiative cooling sets in relatively 
early, and radiative loss can be a significant part of the initial thermal and kinetic energy. Conversely, if 
tc/tf is initially large, the time for the radiative loss is long, and when the gas finally cools radiatively, much 
of its internal energy has been used to overcome gravitation, so radiative loss is not a large fraction of the 
initial internal energy. For figures la and lb, the energy lost to radiation is about 8% and 11% of the total 
internal energy in equation (2.19), respectively. 

When the gas begins losing energy radiatively, it does so very rapidly owing to the increase of the 
cooling rate with decreasing temperature. The pressure of the gas drops rapidly, and now gravitation 
becomes dominant. At this point, because the entropy decreases outwardly as the gas cools radiatively, 
convective and thermal instabilities (Field 1965) ensue, and cloud formation follows. The criterion and 
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some general properties of the convective and thermal instabilities in cooling flows have been worked out by 
Balbus & Soker (1989), but applications to the present problem require more careful and detailed analysis 
since supersonic motions may overrun the instabilities. We shall leave to a future paper the detailed stability 
analyses of the supersonic flow when radiative cooling becomes dominant. For the present purpose, however, 
it suffices to assume that, once radiative cooling sets in and the temperature of the gas drops precipitously 
(e.g., below 10® K) (resulting from radiative cooling) cloud formation ensues, and it completes at cooling 
radius Xc where the temperature drops to ~ lO'*. Note that radiative cooling is a prerequisite for cloud 
formation here; gas can still cool to a very low temperature via adiabatic expansion alone, but in this case 
the flow is stable and the resultant flow is a wind, similar to stellar winds. The onset of radiative cooling is 
signalized by a precipitous drop oitc/tj below unity. 



1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 0.5 
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log X log X 



Fig. 1c 

Figure 1. a). Clockwise from top left, Mach number, dimensionless temperature, normalized 
velocity, normalized pressure, normalized density, and ratio of the cooling time to the Bow 
time, as a function of the dimensionless radius. The initial conditions are Xi = 10"'* and 
Wi = 25. b). Same as in a), but with Xi = 0.6 and Wi = 1. c). Same as in a), but with 
Xi = 10 and Wi = I. 

One may associate the radius at which the temperature drops precipitously because of radiative loss with 
the location of cloud formation. Once clouds are formed, they drop out of the flow and move ballistically. 
The newly formed clouds retain the kinetic energy of the gas at the point of formation. Because the gas loses 
energy to radiation within a time much shorter than the flow time, it has moved little during the period of 
the rapid radiative cooling. Consequently, the velocity, which is mainly affected by the potential and thus 
the relative position of the gas in the halo, remains almost unchanged at the point of cloud formation. We 
may adopt the flnal velocity of the gas (where radiative cooling becomes dominant) as the initial velocity of 
the newly formed clouds. Whether these cooled clouds will eventually escape the galaxy or fall back under 
the influence of the galactic potential depends on the initial cloud velocity relative to the escape velocity, 
which in turn depends on the galaxy mass distribution that gives rise to the rotation curve. 

If the mass distribution (and thus the flat rotation curve) is truncated at radius Tmax, the resultant 
gravitational potential is given by equation (A2). In terms of the dimensionless radius, the escape velocity 
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then 

I 1/2 



V2t;o + • (2.20) 



For example, in figures la and lb the radii at which the temperature of the gas drops precipitously because 
of the radiative loss are Xc = 5.5 x 10"'* and Xc = 2.8, and clouds formed at these radii have velocities 9.4 
and 0.68 times vq, respectively. If x^ax = ^Qxi (or equivalently Vmax = 50rj), the velocities correspond to 
about 370% and 26% of the escape velocity, respectively, so the clouds will be confined in the halo in figure 
lb but will leave the galaxy in figure la. If the gas cannot escape from the galaxy, the final radius to which 
the clouds can coast before being pulled back by gravitation can be estimated from 



Xf = Xc exp 

if X f ^ '^'ma.ri or 



1 (Vc^ 

2 Uo 



(2.21) 



f — o / \ o '^rnaxi (2.22) 

'\ixj> Xmax, where Vc is the initial velocity of the clouds at cooling radius Xc- In figure lb, we have Xf = 3.5, 
or about 30% more distant than Xc- 

Also plotted in figures la-lc are the density, in units of po = A/(47rrgt;o), and the pressure, in units 
of po = vqX/ [AttvI) . The density, after a sharp initial drop as a result of the velocity increase, decreases 
very roughly as x~^ during adiabatic expansion since the velocity then does not vary much. In figure Ic, 
where the gas cannot cool radiatively, the density at the end of the outfiow has a small jump as the fiow is 
stopped by gravitation. A similar jump is also observed in the pressure. In contrast, in figures la and lb, the 
pressure decreases sharply at the final stage of the outfiow because the gas temperature drops precipitously 
as a result of radiative cooling. 

So far we have assumed that the radiative cooling occurs within the edge of the mass distribution of 
the galaxy, i.e., Xc < Xmax- If Xc > x^ax, the gravitational force suffers a jump during the outfiow; inside 
Xmax gravitational force is given by the logarithmic potential in equation (A2), and outside Vmax given by 
the potential in equation (Al) with n = 1. In this case, one can integrate fiow equations (2.14) and (2.15) 
up to Xmax, as done above. Switching to the n = 1 potential at rmax, and adopting as the initial conditions 
the values of the first part solutions at x^ax, we can continue integration of the same equations from Vmax 
outward. (A careful matching at x^ax in terms of the dimensionless variables x and w is necessary because 
the length and velocity scales in eq.[2.10] and [2.11] are different for different values of n.) The effects of mass 
distribution of galaxies can then be investigated by varying x^ax- The calculations are straightforward, and 
for clarity we only summarize the main results here. 

If we denote Xc by the cooling radius obtained by assuming Xc < x^ax, we find that Xc for a given Wi 
is hardly affected by decreasing x^ax below Xc- This is expected because whether the gas can radiatively 
cool depends on the ratio of the radiative cooling time and the fiow time, equation (2.18). The former is 
primarily determined by the initial conditions, and the latter is not much affected by gravitation during the 
adiabatic expansion in a supersonic fiow. Thus for example, in figures la and lb, Xc is increased only by 
0.4% and 27%, respectively, even if we redu ce Xyyiq^x '^ij i-e.j no jnGsswe dovh HgIos for QGloxies . Changing 
^max, however, does affect the fate of the cooled gas since it determines the escape velocity and, hence, 
whether the cooled gas will stay in or leave the galaxy. Similarly, it also infiuences the nature of the outfiow 
if radiative cooling cannot occur, because Vmax in effect determines the total energy (including the potential 
energy) of the gas in the fiow. As seen in the analytical solutions for outfiows without radiative cooling in 
the Appendix, the topology of the fiows is specified by the sign of the total energy. For example, in figure Ic, 
where radiative cooling is not important, if x^ax is reduced below about x^ax — 27, the outfiow topology is 
changed; in this case the total energy of the fiow becomes positive, and the gas leaves the galaxy. The Mach 
number now behaves as that in figure A2 in the Appendix. In the following, we shall be mainly concerned 
with the gas that can cool radiatively; the outfiow solutions without radiative cooling can be readily found 
in the Appendix. We ignore the slight change of the cooling radius Xc caused by varying rmax- However, we 
take into account the effect of finite Vmax on how far the cooled clouds can coast in the halo using equations 
(2.21)-(2.22). 

2d.) General Properties of the Outflow 

From the previous examples, we see that once radiative cooling sets in, the gas temperature drops 
precipitously almost at a fixed radius. This cooling radius Xc for given initial dimensionless temperature w'f 
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is shown in figure 2 for several values of the dimensionless radius a;, = 10 ®, 10 10 ^,10 ^,10 ^, and 1 
(solid curves). (Numerically, Xc is determined when tc/tf in equation [2.18] drops below 0.1.) In the limit 
wf 0, the gas is trapped in the host galaxy, and Xc is equal to Xi. As wf increases, Xc increases, and 
asymptotically approaches cx wf ^^^ for large wf . The transition from the gradual to the asymptotic rise of 
Xc occurs at initial tc/tf ~ 1; below values of wf that correspond to tc/tf ~ 1 at Xi, the gas cools as soon 
as it leaves the base of the galaxy, whereas above those values the gas streams out adiabatically first before 
cooling radiatively. The asymptotic behavior at large wf can be understood as the following. Before the 
onset of radiative cooling the outfiow is approximately adiabatic. From equations (2.14) and (2.15), we find 
approximately 

M2 ~ (x/xi)^/^; (2.23) 
wf{x/x,)-^/', (2.24) 



if radiative cooling and gravity (the first and third terms in the numerators in eq. [2.14] and [2.15]) can be 
ignored. Radiative cooling becomes dominant when tc/tf ^ 1. So substituting equations (2.23) and (2.24) 
into equation (2.18), we find that the cooling radius asymptotically approaches 

4(l-g) 6(2-g) 

XcO,x,'-'- w,'-'- =xj^^^wf ''\ (2.25) 

where the last equality follows from q = —0.6. 

As seen in figures la and lb, even if the radiative cooling time is longer than the fiow time initially, 
adiabatic expansion brings the gas temperature down, and with the assumed cooling function increasing 
with decreasing temperatures, radiative cooling takes over eventually. The gas then cools rapidly unless it is 
stopped prematurely by gravitation (e.g., figure Ic). In practice, however, if the temperature drop through 
adiabatic expansion is too large, our assumed cooling function may not be applicable. For example, if the 
gas cools adiabatically below lO^K before radiative cooling can set in, the rate of radiative cooling for such 
a gas is smaller than the simple extrapolation from the higher temperatures given by equation (2.5). In 
fact, q becomes positive (~ 0.5 between 10® and lO'* K) so the radiative cooling time is increased as the 
temperature decreases. As a result, the gas cannot cool radiatively. In this case, if the initial temperature 
of the gas is above the escape temperature, then the outfiow becomes a wind, a situation similar to stellar 
winds. On the other hand, if the initial temperature is below the escape temperature, the gas remains in 
the galaxy, resulting in a hot corona. 

The cooling function in equation (2.5) is roughly applicable within 10® — 4 x lO^K. The initial gas 
temperature of interest is above a few times 10® K based on X-ray observations (see §3). For example, if 
the initial temperature is 10^ K, and radiative cooling cannot dominate before the adiabatic cooling brings 
the gas temperature down by a factor of 100, the gas will never cool radiatively. From equation (2.24), 
we find that the corresponding upper limits on the calculated cooling radius is about a factor of 32 times 
the initial radius; the cooling radius obtained beyond this limit is fortuitous because the adopted cooling 
function cannot apply here. In figure 2, we plot the values of wf above which radiative cooling cannot become 
dominant before the gas temperature drops by a factor of 100 through adiabatic expansion (shaded area 
above the dash-dotted line). If the initial temperature is higher (lower) than 10^ K, the corresponding w'f for 
given Xi is larger (smaller) than shown in figure 2. If the initial temperature is a x 10^ K, the corresponding 
limit on wf changes by a factor of Q;(i~'*9)/4(2-g) _ q,o 33 compared with the 10^ K case from equations (2.24) 
and (2.25). The above upper limits on wf can also be understood in terms of tc/tf; the limits correspond to 
tc/tf = 42 at the base of the outfiow. That is, if the initial tc/tf is larger, the radiative cooling time is so 
long that radiative cooling cannot become important. 

The region of initial conditions [xi, wf) that results in wind is denoted by WINDS in figure 2; it is 
bounded below by the requirement that wf be above the escape temperature, so the gas can escape the 
galaxy. In the absence of radiative cooling the initial temperature necessary for the gas to escape is 

where /j, is the mean molecular weight, is the mass of hydrogen, and (l){xi) is the potential at Xi. If we 
assume x^ax = ^Oxi, the wind region should be above wf = wf^^ = 2.456 (dashed horizontal line) in figure 
2. Below the wind region the gas cannot leave the galaxy, resulting in a hot corona. 
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Figure 2. Initial dimensionless temperatures versus the dimensionless radius. The solid 
curves represent the cooling radius, Xc, for given wf , calculated for several values of the 
initial radius Xi, from left to right: Xi = 10~®, 10~^, 10~^, 10~^, 10~^ and 1. The 
dotted lines denote the final radius Xf which clouds (formed at Xc) can reach for given 
wf . The initial radii are the same as for the solid curves. Above certain values of wf , Xf 
becomes infinite, indicating that the clouds escape the galaxies. For all the cases, the mass 
distribution of the galaxies is truncated at x^ax = ^Qxi. The dash-dotted line shows the 
boundary of wf as a function of Xi above which the initial conditions are such that radiative 
cooling cannot be important throughout the Bow (shaded area). The shaded area is further 
divided into two regions; above the escape temperature wf = 2.456 (dashed horizontal line), 
the outflow results in galactic winds (denoted WINDS), and below it the outEows lead to 
coronae. The other shaded area in the lower right, denoted GRAVITY, represents the region 
of initial condition (xi, wf) where the gravitational potential of the galaxies dominates and 
the gas is stopped by gravity before it has a chance to cool radiatively; in this case, the 
outEows also result in galactic coronae. The two crosses show the range of initial conditions 
(xi, wf) for dwarf galaxies with rotation velocity Vdr = 50 km/s (left cross) and normal 
galaxies with rotation velocity Vdr = 225 km/s (right cross). The vertical extension of 
the crosses represents the initial temperature range considered, Ti = iQ^^^^o.s^ j-^g 
horizontal stretch of the crosses represents the range ri kpc/^i = 10^^^ for both dwarf and 
normal galaxies. 

For Xi ^ 1, gravitational force is strong relative to radiative cooling as discussed in §2a, so the gas will 
be pulled back to the galaxy before it has a chance to cool radiatively (cf., figure Ic). In figure 2, we show the 
region of (xi, w'f) where the gas is confined to the halo by gravitation and radiative cooling is not important 
during the fiow time. The region is bounded below by tc/tj ~ 1 since obviously radiative cooling must not 
be too strong for gas to rise at all. It is also bounded above by the escape temperature, as the gas can leave 
the galaxy in this case. The escape velocity is taken to be that given by x^ax = ^Qxi. 

We have seen in §2c that in case the gas cools radiatively and the temperature drops precipitously at the 
end, the cooled gas can coast farther beyond the cooling radius Xc because the newly formed clouds inherit 
the finite kinetic energy from the hot gas. In figure 2 we show for given w'f the maximum radius, xj, which 
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the cold clouds can eventually reach, based on equations (2.21) and (2.22) (dotted curves). Calculations 
are made for several initial radii Xi = 10~®, 10"'*, 10~^, 10~^, 10~^, and 1. For definiteness, we take 
Xmax = ^Qxi] in general, if < x^ax, Xf is independent of x^ax, and ifxf > x^ax, Xf is proportional to 
Xmax, (eqs. [2.21] and [2.22]). An indication of the relative importance of radiative loss can be found from 
the minimum wf that allows the cooled gas to escape the galaxy. For example, in figure 2 we find that for 
Xi = 10~® the initial temperature required for escape is wf ^ 3.9. The dimensionless escape temperature 
from equation (2.26) is Wg^g = 2.456 for x^ax = ^Oxi. Thus in case Wi = 3.9, about 60% of the initial 
internal energy is lost to radiation. As Xi increases, the relative magnitude of radiative cooling diminishes. 
Much of the energy is consumed in overcoming gravitation, so radiative loss is only a small fraction of the 
initial internal energy. For instance, from figure 2 we note that for Xi = 0.1 the required wf for the escape 
is ^ 2.55. In this case, we infer that only about 4% of the initial internal energy is lost to radiation. 

3. APPLICATIONS 

In this section, we switch to dimensional variables and make attempts to connect with various obser- 
vations. More specifically, we first examine the difference in characteristics of the cooling outfiows between 
dwarf and normal galaxies. We then calculate the X-ray surface brightness, the mean gas temperature, and 
the mean surface brightness of the O VI emission lines for some representative outfiow solutions. 

3a). Dwarf and Normal Galaxies 

The cooling function for a gas of nearly cosmic abundance at 10^ < Tg < 4 x 10^ K may be taken 
approximately as (e.g., Sutherland & Dopita 1993) 

q = -0.6 

A = 3.4 X lO^^'cgs 

For gas with lower abundance, the above q may still be roughly applied, but A is smaller. For example, A 
is a factor of about 3 (5) smaller for a gas of 1/10 (1/100) the cosmic abundance (e.g., Bohringer & Hensler 
1989; Sutherland & Dopita 1993). Non equilibrium of collisional ionization may also reduce the cooling rate 
at a temperature below about 4 x 10® K. The reduction may be up to a factor of 4 at 10® K (e.g., Shapiro 
& Moore 1976; Schmutzler & Tscharnuter 1993; Sutherland & Dopita 1993). If a cooling rate other than 
equation (3.1) is to be used, the corresponding results are modified according to figure 2; x scales with AX 
(so reducing A is equivalent to increasing A). 

Assuming a fiat rotation curve (n = 0), from equations (2.10a) and (2.11a), we obtain 

Ai 

ro = 278-gYkpc, 

^100 

where Ai = X/ [Mq /yY), and t;ioo = i'o/(100km/s). Thus the dimensionless radius is 

^ = 3.6 X 10-3^^%^, 
Ai 

and the dimensionless temperature is 

^.2 = 2.2-^, (3.4) 

"lOO 

where rt-pc is the radius in kpc, and Tq = Tg/10® K. Note that equation (3.3) only depends on ratio rt-pd^i- 
Now apply equations (3.3)-(3.4) to initial conditions Xi and w'f where r = r,- and Tg = T,-. The initial gas 
temperature is primarily determined by the heating of the ISM by supernova remnants in the galaxy, which 
depends on the internal star formation rate per unit volume (e.g., McKee 1990). Observationally, the star 
formation rate is found to correlate more with morphological types rather than the mass of galaxies (e.g., 
Sandage 1986), so one may expect that the initial temperature is not a strong function of the galaxy mass. 
However, Xi is a strong function of the circular velocity, and, thus of the mass of the galaxy. As a result, for 
a given mass-loss rate, massive galaxies in general have large Xi and small Wi, and dwarf galaxies in general 
cluster around small Xi and large Wi. 

As examples, in figure 2 we show the inferred [xi, w'f) for dwarf galaxies with a circular velocity of 50 
km/s if Ti_6 = 3 and ri^kpc/^i = 1 (e.g., n^kpc = 0.3 and Ai = 0.3), where Ti_6 = T^/IO^K and ri^kpc = n/kpc. 
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(3.1) 



(3.2) 



(3.3) 



The cross on the left in figure 2 represents the uncertainties of the various parameters; the vertical extension 
of the cross results from the adopted temperature range Ti e = 10'^^®^'^'®, and the horizontal stretch refiects 
the adopted of ri kpc/^i = W^^^^ (e-g-, Ai = 0.01 — 10 for fixed ri kpc = 0.3). Similarly, we also present the 
initial conditions expected for normal galaxies like our own with circular velocity v^r = 225 km/s if Tj_6 = 3, 
ri,kpc/^i = 1 (e.g., ri kpc = 3 and Ai = 3). A similar range of the parameters is taken: Ti e = lO*^'®^*^'®, and 
ri,kpc/^i = 10=*=^^® (e.g., Ai = 0.1 - 100 for fixed ri^kpc = 3) (right cross). 

A comparison with the previous calculations of the cooling radius Xc and the final radius Xf for the 
cooled clouds in figure 2 shows that dwarf galaxies in general should have halos of hot gas many times their 
initial radii. Depending on the initial temperature and density, the gas can either cool or result in galactic 
winds. In either case, the gas is most likely to leave the galaxy because the temperature is usually above 
the escape temperature. In contrast, for massive galaxies like our own, unless the temperature is extremely 
high, the gas is in general confined within the galaxy. If the gas density is high enough (large A and thus 
small Xi), radiative cooling becomes important at the end and the gas cools in the halo. But typically the 
cooled gas is unable to travel much farther from the cooling radius because of the strong potential. If the 
gas temperature is below the escape temperature, and the density is low, the gas will be pulled back to the 
galaxy before it has a chance to cool radiatively. In this case, the fiow will result in a galactic corona. 

3b). Surface Brightness and Temperature Structure 

We show in figure 3 some examples of the temperature structures and density profiles in the cooling 
outfiows from dwarf galaxies with Vdr = 50 km/s, and from normal galaxies with Vdr = 225 km/s. We 
take ri kpc roughly as the region of active star formation in galaxies. We consider a modest outfiow with 
Tj- = 3 X 10® K, and a strong outfiow with T,- = 2 x 10^ K. The input parameters and the general results of 
the outfiow calculations are summarized in Table 1. 

We also present the surface brightness in 0.1-2.2 keV (roughly the ROSAT energy range) and in 1.6-8.3 
keV bands, and the total surface brightness. The 1.6-8.3 keV hard band may be compared with the X-ray 
observations by GINGA (2-20 keV), BBXRT(0.5-10 keV) and AS CA (0.5-10 keV). The surface brightness is 
calculated through the integral T{r) = J^°° z p'^ e[Tg) / z'-^ — r'^)dz, where the emissivity of the gas, e{Tg), 
is taken from Suchkov et al. (1994). 

If the outfiows are dominated by adiabatic expansion, the temperature is proportional to r~^^^, and 
the density is roughly p cx r~^, so the surface brightness is ~ rp'^e cx r3-4q/3 _ j,-2.2^ p^j, ^j^g modest 
outfiows considered above, the initial cooling time is slightly larger than the fiow time {tc/tj = 5), so 
radiative cooling soon becomes important as the fiow streams out. As a result, the total surface brightness is 
much shallower than r~^^ because the temperature drops more rapidly than r~^l^ . The surface brightness 
in the ROSAT band, however, drops much more sharply than the total surface brightness because of the 
temperature decrease. The gas is also too cool to be easily seen in the hard energy band. In contrast, for 
the strong outfiow, the initial cooling time is relatively large compared with the fiow time {tc/tj = 20), so 
the fiow first expands adiabatically, and the total surface brightness decreases roughly as r~^^. At the end 
of the fiow, however, the cooling time becomes much shorter than the fiow time, and the temperature drops 
rapidly as a result of radiative cooling. Consequently, the total surface brightness decreases less steeply than 
in the adiabatic phase before the final precipitous decrease. 

The cooled clouds formed in both the modest and strong outfiows from dwarf galaxies have large initial 
velocities well beyond the escape velocity for any reasonable mass distribution, so they leave the galaxies. 
For the normal galaxies, the escape velocities are 461, 514, 578, and 613 km/s if the mass distribution of the 
galaxies is truncated at 30, 50, 100, and 150 kpc, respectively. From Table 1, we thus see that the cooled gas 
is most likely to stay in the galaxies in the modest outfiows, and it leaves galaxies in the strong outfiows. 

The calculations of the modest outfiow from dwarf galaxies may be compared with the ROSAT ob- 
servations of NGC 1569 {vcir = 50 km/s) (Heckman et al. 1995). The extent of the observed X-ray halo 
emission is about several hundreds pc in radius with a luminosity in the ROSAT band of about 4 x 10^® 
erg/s. The Ha emission in NGC 1569 has an extent of about 200-300 pc in the central region (Walker 1991). 
The calculations of the modest fiows from normal galaxies may be compared with the ROSAT observations 
of NGC 891 (Bregman & Pildis 1994). The full extent of the observed X-ray halo emission is about 7 kpc 
in radius with a luminosity in the ROSAT band of about 4 x 10^® erg/s. The scale height of the extended 
Ha emission in the central region of NGC 891 is about 3-4 kpc (Rand, Kulkarni, & Hester 1990; Dettmar 
1990). Extended X-ray emissions from a few other galaxies have also been detected with ROSAT, including 
M82 (Bregman, Schulman, & Tomisaka 1995) and NGC 4631 (Wang et al. 1995). We shall present more 
detailed comparisons with observations elsewhere. 
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To compare with the average temperature of the entire halo derived from X-ray observations, we cal- 
culate the mean temperature averaged over the outflow region and weighted by the emission: < T >= 
J Airr^CTgdr I [ j Airr^ C dr) . Figure 4 shows < T > as a function of the initial temperature Ti for the four 
examples shown in flgure 3: Not surprisingly, the mean temperature is systematically lower than the initial 
gas temperature (by a factor of 2-6) because of the cooling of the gas away from the galaxies. Correctly 
taking this into account is important in understanding various physical conditions (such as the heating of 
the gas) in the host galaxies. 




r(kpc) r(kpc) 

Fig. 3 



Figure 3 Gas temperature (solid curves), gas density (dotted), total surface brightness 
(dashed), surface brightness in 0.1-2.2 keV (long-dashed), and surface brightness in 1.6- 
8.3 keV (dot-dashed) as a function of the radius for some representative cooling outEows. 
The level of the surface brightness is labeled on the right axes. a). Results of the outflow 
from dwarf galaxies (vdr = 50 km/s), with initial temperature Te = 3 and mass-loss rate 
Ai = 0.3. h). Same as a), but with Ti = 20 and Ai = 10. c). Results of the outflow from 
normal galaxies (vdr = 225 km/s), with Te = 3 and Ai = 3. d). Same as c), but with 
Ti = 20 and Ai = 100. Other model parameters can be found in Table 1. In a) and c), the 
fluxes in 0.1-2.2 keV are below the flux ranges shown. 

The cooling gas can also be detected via emission lines, such as the O VI doublets (A1034 A). For 
examples, for a gas cooling from Te = 3 (Tg = 20) down to Tg = 0.1, the energy emitted in the VI 
emission lines is about 3% (0.5%) of the total energy lost to radiation (Voit, Donahue, & Slavin 1994). Table 
1 gives the mean surface brightness in O VI, calculated as the surface brightness averaged over the flow 
region inside re] the actual surface brightness in the inner regions should be much higher than the average. 
The predicted mean surface brightness may be within the reach of current UV experiments. For example, 
the detectable O VI line flux (3 cr) of the spectrometer recently flown on a rocket is 2 x 10~^ erg/s/cm^/sr 
(Rasmussen & Martin 1992). The Hopkins Ultraviolet Telescope (HUT) flown aboard the Space Shuttle on 
the Astro-1 mission also has a similar sensitivity for a typical exposure time (~ 2000 s) (Davidsen et al. 
1992). Since HUT on Astro-2 mission scheduled next year is expected to improve the sensitivity by a factor 
of few (W. P. Blair 1994, private communication), the outflows, especially from dwarf galaxies, should be 
detectable by HUT. 
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Figure 4. Mean temperature averaged over the entire Bow region as a function of the initial 
temperature (eq. [3.7]), for the outEows shown in figures 3a (solid curve), 3b (dotted), 3c 
(dashed) and 3d (dot-dashed) . 

4. DISCUSSION AND CONCLUSION 

We study the outflows of hot gas from galaxies. In particular, we investigate the general properties 
of gas outflows. Emphasis is made on understanding the outflow topology, the extent of the outflow, the 
efficiency of radiative cooling, and the fate of the cooled gas. As a flrst step, we consider problems of 
steady, radial outflows with initial conditions specifled at the base where the flow originates. Gravitational 
potentials in galaxies are assumed to have either power-law or logarithmic dependence on the radius. We 
demonstrate in analogy with stellar wind problems that without efficient heating the outflow cannot overcome 
the sonic barrier, so the outflow remains either subsonic or supersonic (Appendix). While subsonic flows 
may be subject to various convective and thermal instabilities and may not travel much farther from the 
base because of the relatively small velocities, the branch of supersonic solutions provides the most promising 
route for gas to reach large distances away from the galaxies. 

We approximate the cooling function by a power-law of the gas temperature. Along with the assumed 
forms of the galactic potential, there exist two natural scales in the problem (eq. [2.10] and [2.11]), and the 
flow equations can be reduced to dimensionless form (eq. [2.14] and [2.15]). The outflow problems then can 
be formulated in terms of solving ordinary differential equations with two-parameter initial conditions: the 
dimensionless initial radius and the dimensionless initial temperature. This greatly reduces the parameter 
space one would need to explore. The dimensionless radius x (eq. [2.12]) contains the information concerning 
the importance of the gravitational potential relative to radiative cooling, and the dimensionless temperature 
(eq. [2.13]) measures the thermal energy of the gas relative to the potential energy. 

We obtain supersonic solutions for the various initial conditions representing a variety of galaxies. A few 
examples of those solutions are shown in flgure 1. The characteristics of the outflows are determined by the 
ratio of the cooling time to the flow time, tc/tf (eq. [2.18]), and by the initial dimensionless radius Xi. If the 
galactic gravitational potential is not important (i.e., Xi ^ 1), the outflows are distinguished by the initial 
values oitc/tj: (1) If initially tc/^/ ~ 1, the gas cools as soon as it streams out of the galaxies. (2) If initially 
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tc/i} ~ 1, gas first cools adiabatically, and radiative cooling becomes dominant later on. (3) Radiative 
cooling cannot become important iftc/tf ^ 1 (e.g., tc/tf ^ 40 for a gas with an initial temperature 10^ K) 
because the cooling time becomes too long. In the last case, the outflows either become galactic winds or 
stay in the galaxy as a corona, depending on whether the initial temperature is above the escape temperature 
of the galaxy. If, however, gravity is strong (i.e., Xi ^ 1), the gas is pulled back to the galaxy before it has 
an opportunity to cool radiatively, and the outflow results in a hot corona. It is important to note that 
tc/tf ~ 1 at the base radius may not be used as the sole criterion for the onset of the radiative cooling; 
adiabatic expansion brings the gas temperature down, so tc/tf decreases in the flow. Only when tc/tf 1 
can radiative cooling be ignored. 

In case radiative cooling does become dominant, the gas temperature drops precipitously owing to the 
increase of the cooling rate with decreasing temperature. One may define a cooling radius where the gas 
temperature drops rapidly because of radiative cooling. Once radiative cooling sets in, the gas becomes 
unstable to convective and thermal instabilities, and the cloud formation ensues. The newly formed clouds 
inherit the kinetic energy of the gas at the point of cloud formation and coast farther from the cooling radius. 
For large initial temperature and shallow potential wells, the clouds escape the galaxy, while for galaxies 
with deep potentials, the clouds remain in the halo. Our general results on the outfiows and the fate of the 
cooled gas are summarized in figure 2. 

We apply our results to dwarf galaxies and normal galaxies. For reasonable ranges of the parameters, we 
find that dwarf galaxies generally have extended X-ray halos, and the outfiows leave the galaxies eventually 
because of the relatively shallow potential wells. The gas may cool and forms clouds, in which case the 
cold clouds also leave the galaxy. If the gas density is low, radiative cooling never becomes important, in 
which case the outfiow is in the form of galactic winds. In contrast, for large galaxies like our own, gravity 
is relatively strong, and the gas may not have a chance to cool before being pulled back to the galaxy. Even 
if the gas density is high enough so that radiative cooling becomes important eventually, the newly formed 
clouds are most likely to remain in the galaxy because of the strong gravitational field. In the extreme case 
where the temperature is larger than the escape temperature and radiative cooling is not important, galactic 
winds take place. We summarize the likely initial conditions for dwarf and normal galaxies in figure 2. 

Some representative radial profiles of the gas temperature, the density, and the surface brightness of the 
cooling outfiows are presented for dwarf and normal galaxies (figure 3). We find that the extent of the X-ray 
emission in the halo is usually much smaller than the outfiow regions. The total radiation from the outfiow 
that can be detected in the ROSAT band is about one order of magnitude smaller than the total radiation 
available for the initial temperature of about 3 x 10® K, and the fraction of the radiation emitted in the 
ROSAT band increases to about 30% of the total radiation if the initial temperature increases to 2 x 10^ K. 
The X-ray in the hard (e.g., ASCA) energy band is smaller still, and we expect only the high temperature 
outfiows [Ti ^ 10^ K) may be detected in the hard band. We calculate the mean temperature averaged over 
the entire outfiow region, and we show that the mean temperature is usually a factor of 2-6 lower than the 
base temperature owing to the cooling of the gas away from the base (figure 4). This difference is important 
in understanding the physical conditions in the host galaxies, as usually only the mean gas temperature of 
the entire halo can be derived from observations. We also estimate the strength of the O VI emission lines 
(A = 1034 A). The mean surface brightness of the O VI doublets, especially from outfiows in dwarf galaxies, 
is within the sensitivity of current UV experiments. We summarize the results of the calculations for the 
representative outfiows in Table 1. 

As mentioned in §1, the steady-outfiow calculations may be applied to galaxies where the heating lasts 
longer than the fiow time. However, for some extremely strong starburst galaxies, the fiow region can be 
so large that the finite duration of the heating must be considered. An example is NGC 3628 (vdr = 240 
km/s) which has been detected by ROSAT to have an extended X-ray halo with radius of ~ 30 kpc (Dahlem 
et al. 1995). Since the galaxy may have experienced a starburst in the recent past (Fabbiano et al. 1990), 
the last period of heating by the burst of supernovae may be shorter than the crossing time for the outfiow 
region. If so, time-dependent calculations are necessary. Also, since we have assumed spherical symmetry, 
our calculations cannot model nonradial outfiows. 

Our results may have significant implications for observations of QSO absorption line systems. More 
detailed discussion can be found in Wang (1995); here we only briefiy mention possible applications of 
the model calculations presented above. Recent observations of the absorption lines in QSO spectra have 
increasingly shown evidence for the absorbing gas at large distances away from the associated galaxies, in 
particular for the heavy element absorbing gas, which has long been suspected to reside in galactic halos 
(Bahcall & Spitzer 1969). The identification of the galaxies causing the Mg II absorption lines in the 
intermediate redshifts show that the absorbing gas has a radius around the associated galaxies of about 50 
kpc (e.g., Bergeron & Boisse 1991; Yanny et al. 1990; Bechtold & Ellingson 1992; Steidel 1993). Since the 
absorption lines of higher ionization states such as C IV are even more numerous (e.g., Sargent, Boksenberg, 
& Steidel 1988), this implies an even larger size for their absorbing regions if the gas is associated with 
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galaxies. More recently, studies of the relation of the low column Lyman a forests with nearby galaxies 
show that they may be spatially correlated with the galaxies (Morris et al. 1993; Lanzetta et al. 1995), 
indicating a possible connection between galaxies and the absorption lines. If indeed Lyman a forests are 
also associated with galaxies, the size of the absorbing gas is even larger than that inferred for the heavy 
element absorbing gas, ^ 160 kpc (Lanzetta et al. 1995). This size in fact agrees with the results of studies 
of Lya lines in pairs of quasars at z = 2. 

The huge size of the absorbing gas associated with the galaxies is rather hard to understand in the 
context of the normal galaxies like our own (Wang 1993). As seen in NGC 891 and from many other failed 
attempts to detect the halo X-ray emissions from normal galaxies, large halos are rare phenomena rather 
than the norm for normal galaxies. Indeed, estimates of the gas temperature of the hot phase in the Galactic 
ISM show that the temperature is T,- ^ 10® K (McKee & Ostriker 1978), and from figure 2 we see that the 
gas is trapped in the halo within a few kpc (also see Houck & Bregman 1990). Less massive galaxies, in 
particular dwarf galaxies, may, however, provide the needed gas at large distances away from galaxies. The 
gas temperature of the hot phase in less massive galaxies may be comparable to that of normal galaxies 
because the same heating mechanism through the shocks of supernova remnants is at work. But their 
gravitational potentials are much shallower than those of normal galaxies, therefore the gas is most likely to 
leave the galaxies. Our present work shows that the outfiows from dwarf galaxies can cool radiatively for 
reasonable mass-loss rates (e.g., a few tenths of MQ/yr). The cooled gas is most likely to leave the galaxies 
as it inherits the kinetic energy of the hot gas which may not be exhausted by the cooling. The cooled gas 
thus can cause absorption over an area much larger than the original size of the galaxies. Note that radiative 
cooling of the gas may be crucial to accounting for both the observed low ionization states of the heavy 
elements and the narrow Lya absorption lines. 

I am especially grateful to Tim Heckman for useful discussions on observations of the gas outfiows 
from galaxies throughout the course of the project. I also thank Mark Voit for helpful suggestions on the 
presentation of the present work. The author also thanks Steve Balbus, Mike Dahlem, and Bill Blair for 
useful discussions. Support from Alan C. Davis Fellowship at the Johns Hopkins University and Space 
Science Telescope Institute is gratefully acknowledged. 
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APPENDIX 

STEADY, SPHERICAL FLOWS WITH CONSTANT 7 



Here we consider the flow problem given by equations (2.1)-(2.3) but without radiative cooling (A = 
or g — — cxd). Unlike in the rest of the paper, the adiabatic index is taken as a constant that can be difli"erent 
from 5/3, emulating possible heating and cooling processes. Gravitational potential is assumed according to 
equation (2.9): 

Hr) = -:^- (n^O) (Al) 



For n = 0, we assume a cut-off radius for the mass distribution in the galaxy, 
divergence, so within rmax the potential is. 



to avoid the logarithmic 



(r) 



In 



(n = 0) 



{A2) 



where v^r is the circular (or rotation) velocity of the galaxy. Beyond rmax the potential is the same as that 
of a point mass and is proportional to 1/r (n = 1 in eq. [Al]). Behaviors of the flow outside rmax are readily 
understood from properties of the flow in a potential given by equation (Al) with n = 1, so for n = we 
shall concentrate on studying the flow properties within Tmax- 

Flow solutions for spherical accretion problems with constant 7 in case n = 1 have been found by Bondi 
(1952), and the corresponding wind solutions have been studied by Parker (1958). A good review can be 
found in Holzer & Axford (1970). Here we extend flow solutions to arbitrary values of n to take into account 
galactic gravitational potentials. 

Upon taking q —00 in equations (2.6) and (2.7), we flnd that M'^ = 1 is a critical point, same as in 
the case with radiative cooling. To proceed smoothly through the sonic point, we require 



2c' - r-^ = 0. 
dr 



(at M' = 1) 



The total energy per unit mass, E, is given by 



E 



7 p 



2 7-lp 

Substituting equation (A3) into (A4), for n we obtain 



+ d>(r) 



{A4) 



E 



B 

nr" 



(4 + n) - (4 - n)7 
4(7-1) 



(A5) 



where the subscript s denotes the critical point where M' = 1 (sonic point). Namely, at the sonic point, the 
total energy is a critical fraction of the potential energy. To pass the sonic point smoothly, we require this 
fraction to be positive (£" > 0). This is true only if 



7 < 



n 5 

- < TT, 

n 3 



(A6) 



the last inequality is valid if n < 1. If n = 1, we obtain 7 < 5/3, recovering the condition for transonic 
flows found by Bondi (1952) in accretion problems and by Parker (1958) in stellar winds. That is, for the 
point-mass potential, an adiabatic flow is either subsonic or supersonic. However, the condition for transonic 
flows is in general more stringent for potentials with exponent n different from unity as readily seen in 
equation (A6). In fact, upon taking n — 0, we flnd that the requirement becomes 7 < 1. (The same result 
can be obtained by using equation [A2] directly). Obviously, flows cannot pass the sonic point smoothly 
even for an isothermal gas in a potential that gives rise to a flat rotation curve. The gravitational held in 
this case is much stronger than that of a point mass, and the nozzle effect necessary for a transonic flow is 
much more difficult to achieve. Thus, to overcome strong gravity, the gas has to be heated with increasing 
distance away from the origin in order to pass the sonic point. 
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X=r/r, 



3 



Fig. A1 

Figure Al. Mach number M as a function of the radius r normalized to the sonic radius 
(eq. [AlO]). We take 7 = 1.25. The solid curves are for the point-mass potential (n = I) 
and the dotted curves are for n = 0.5 potential in equation (Al). The transonic solutions 
passing through the sonic point have S = I (eq. [A9]). The two supersonic (or subsonic) 
solutions for each n have S = 0.99 and 0.95 away from the sonic point. The two solutions 
for each n at x < 1 (or X > IJ that proceed from supersonic to subsonic or vice versa but 
do not pass the sonic point have S = 1.01 and 1.05 away from the sonic point (one of curves 
for n = 0.5 is out of the range to the right). 

To obtain flow solutions, we note that without radiative cooling, equations (2.1)-(2.3) can be integrated 
directly. Following Holzer & Axford (1970), we choose the mass-loss rate. A, the energy per unit mass, E, 
and the speciflc entropy, s = Cy \n[Pl p^) [cy is the speciflc heat capacity at a constant volume), as constants 
of integration.^ The reference quantities are taken at the sonic point. From constant entropy, we may write 
the gas density 

[Al) 



p = ^1 (-'-!) c(-'-i) exp 



(7-1)' 



From the conservation of mass and using equation (A7), we obtain the sound speed in terms of the sound 
speed at the sonic point: 

2 

_ 2(t-1) _ 4(t-1) 

= SM (-'+!) X , (-48) 



where x = ^'l^'s, and 



5 = exp 



2(s-s.) 



(A9) 



Again, subscript s denotes the variables at the sonic point. Substituting equation (A9) into equations (A4) 



An alternative approach is to use Bernoulli's equation directly (Bondi 1952). 
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and using equation (A5), we obtain solutions for transonic flows in terms of the Mach number for n ^ 0: 



(4 + n) - (4 - n)7 
4(7-1) 



+ X~ 



'(T-i) , / M2 1 
V (t+1) = S 1 

^ * 2 7 - 1 



2(t-1) 

M (-'+1) 



{AW) 



Note that solutions for 7 = 1 can be obtained by taking the (7 — 1) root on both sides of equation (AlO); 
letting 7—^1 and using the definition (1 + rjxY^'^ = in the limit ?y — 0, we have 



■Alnx = - InM^ + C\, 



(All) 



where C'l is a constant of integration. 




Fig. A2 

Figure A2. Mach number M for solutions with £" > as a function of the radius r normalized 
to the chosen reference radius ri (see the text). The solid curves represent the Bow solutions 
for the n = 1 potential in equation (Al) (the point-mass potential) . We take 7 = 5/3. The 
values of 5 are from left to right 0.300, 0.425, 0.550, 0.675, and 0.800. The dotted curves 
denote the Bow solutions with E > in isothermal spheres (eq. [A2]), and we have taken 
Xmax = 3. The values of 5 are the same as those for the solid curves. 

For isothermal spheres n = 0, similar transonic solutions can be obtained if the sonic point is inside 
4. 



(I + 7) 
4(7-1) 



Inx 



M2 



7 



1 



2(-,-l) 

M (-'+!) 



(A12) 



Some examples of transonic solutions (AlO) with n = 1 and 0.5 are given in figure Al; transonic solutions 
(A12) for isothermal spheres with 7 < 1 are similar, so we do not plot them separately here. Once the Mach 



■* If the sonic point is outside rmax, the solutions are topologically similar to those in n = 1 potentials 
given by equation (Al) because outside Vmax gravitational field is equivalent to that by a point with a mass 
equal to the total mass within rmax- 



22 



number is known, it is straightforward to obtain the temperature (cx and the density from equations 
(A8) and (A7). 

If equations (A6) cannot be satisfied, transonic solutions are not possible. In this case, we may choose 
some different reference quantities other than at the sonic point, e.g., at r = ri (see below). From equation 
(A4), and using equations (A7) and (A8), we obtain for n ^ 



E 



s 



7 



M (-'+!) 



(AU) 



where Vq = B/nr" . If £" > 0, we choose Vq = E, and if £" < 0, we take Vq 
arbitrarily. Similar solutions can also be obtained for n = 0: 



-E. If £" = 0, Vq can be chosen 



4+1 



In 



X 



Xr. 



X 



(t+1) 



s 



1 



M (-'+1) 



where Xr. 



I'max/i'i, Vq = f cjV ' ^^'^ have assumed ri < r-max- If -E" > 0, we choose E 



(AU) 



and if 



^ < 0, we choose E — — VQlu^fyiQx • 

Beyond rmax (if solutions do exist), one may extend equation (A14) 
with solution (A13) with appropriate boundary conditions matched. 





0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

X=r/rj X=r/ri 

(a) (b) 
Fig. A3 

Figure A3, a). Mach number M for solutions with E < as a function of the radius r 
normalized to the maximum radius ri (vertical lines at x = which the gas can reach. 
Gravitational potential is given by equation (Al) with n = 1 (the point-mass potential). 
We take 7 = 5/3. The values of S are from right to left 0.100, 0.175, 0.250, 0.325, and 0.400. 
h). Same as in a) except that the potential is given by equation (A2) (isothermal spheres). 
The values of 5 from the outer to the inner curves are 0.050, 0.100, 0.140, 0.160, and 0.175. 

Examples of the solutions for either suhsomc or supersonic fiow based on equations (A13) and (A14) 
with E > Q are shown in figure A2. For n = 0, we have taken Xmax = 3, but solutions for other values of 
Xmax are topologically similar. It is clear from equation (A13) that for the particular combination 



(7 + I) 
4(7-1)^ 



(A15) 
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the Mach number is independent of radius for E = 0. For example, an adiabatic gas (7 = 5/3) has a constant 
Mach number in the point-mass potential (n = 1) if £" = 0. Except for the special conditions imposed by 
equation (A15), the behavior of the solutions with E = are similar to those with E > 0. However, solutions 
with E < behave drastically differently from E > solutions. It can be easily inferred from equations 
(A13) and (A14) that ri in this case defines a maximum radius beyond which fiow solutions cannot exist. 
Physically, because the gas has an energy lower than that required to escape, it is trapped in the galaxy. In 
figure A3, we show some examples of solutions (A13) and (A14) for E <0. 

The asymptotic behavior of at x — cxd for both transonic fiows and nontransonic fiows with E > 
can be obtained from equations (AlO) and (A13). In the supersonic regime, we have 

M2^x'*^"'^ (^16) 

and in the subsonic regime we get 

From equations (A7) and (A8), it follows that the density in supersonic fiows asymptotically approaches 
p (X J while both the temperature and the density in subsonic fiows approach a constant. Namely, a 
subsonic fiow requires a confining pressure at infinity. Note that this asymptotic behavior can be also applied 
to isothermal spheres which have to be truncated at Tmax', at r > Tmax, the potential is of n = 1. 
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Table 1 

Model Calculations in §3& 



Dwarf Galaxies Normal Galaxies 



Input 


Vcir (km/s)" 




50 




50 




225 




225 




r; (kpc)* 




0.3 




0.3 




3 




3 




Ti {Ky 


3 


X 10^ 


2 


X 10^ 


3 


X 10^ 


2 


X 10^ 








0.3 




10 




3 




100 


Result 


r-c (kpc)^ 




1.5 




5 




13 




50 




Vc (Vin/s)^ 




490 




1310 




295 




1200 




rj (kpc)^ 




CXD 




CXD 




32 




CXD 




L(total) (erg/s)'' 


6 


X 10^9 


4 


X 1041 


8 


X 1040 


4 


X 1042 




L(0. 1-2.2 keV) (erg/s)' 


4 


X 10^8 


1 


X 1041 


6 


X 10^9 


1 


X 1042 




L(1.6-8.3 keV) (erg/s^ 


2 


X 10^5 


8 


X lO^s 


3 


X 10^6 


8 


X 1040 




L(0 VI) [eYg/sf 


2 


X 10^8 


1 


X 1040 


4 


X 10^8 


3 


X 10^9 




F(0 VI) (erg/s/cmVsr)' 


2 : 


>< 10"^ 


1 


X 10"^ 


7 : 


>< 10-9 


4 


X 10-9 



" Galaxy rotation velocity. 
* The initial (base) radius. 

The initial temperature at the base. 

The mass-loss rate. 

The cooling radius. 

^ The velocity at the cooling radius, and thus the initial velocity of the clouds. 

^ The final radius which the cool clouds can coast to. The mass distribution is truncated at Vmax = 50rj' (15 

kpc and 150 kpc for dwarf and normal galaxies, respectively). 
The total luminosity integrated over the entire outfiow. 
' The luminosity in 0.1-2.2 keV band (corresponding to the ROSAT energy band) integrated over the entire 

outfiow. 

^ The luminosity in 1.6-8.3 keV band (cf., the energy bands of GINGA, BBXRT and ASCA) integrated over 
the entire outfiow. 

^ The total line luminosity from the O VI emission lines integrated over the outfiow region inside Vc- 
' The mean surface brightness of the O VI emission lines average over the outfiow region inside Vc- 
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